Matrix product states for quantum metrology 
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We demonstrate that the optimal states in lossy quantum interferometry may be efficiently sim- 
ulated using low rank matrix product states. We argue that this should be expected in all realistic 

quantum mctrological protocols with uncorrclatcd noise and is related to the elusive nature of the 

Hciscnbcrg precision scaling in presence of decoherence. 
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Over the recent years, advancements in quantum 
engineering have pushed non-classical concepts such 
as entanglement and squeezing, previously regarded 
as largely academic topics, close to practical applica- 
tions. Quantum features of light and atoms helped to 
improve the performance of measuring devices that 
operate in the regime were the precision is limited 
by the fundamental laws of physics [1]. One of the 
most spectacular examples of practical applications 
of quantum metrology can be found in gravitational 
wave detectors [2], were the original idea [3] of em- 
ploying squeezed states of light to improve the sens- 
itivity of an interferometer has found its full scale 
realization [4, 5]. No less impressive are the works of 
the 2012 Nobel Prize winner David Wineland were 
entangled trapped ions have been demonstrated to 
provide spectroscopic resolution enhancement cru- 
cial for the operation of the atomic clocks [6-8]. 

When standard sources of laser light are being 
used, any interferometric experiment may be fully 
described by treating each photon individually and 
claiming that each photon interferes only with itself. 
Sensing a phase delay (f) between the two arms of the 
interferometer via intensity measurements may be 
regarded as many independent repetitions of single 
photon interferometric experiments. N independ- 
ent experiments results in the data that allows the 
parameter (p to be estimated with error scaling as 
1 /VN — the so called standard quantum limit or the 
shot noise limit. If, however, an experiment can- 
not be split into N independent processes, as is 
e.g. the case with the N probing photons being en- 
tangled, the above reasoning is invalid and one can 
in principle achieve the 1/N estimation precision — 
the Heiseberg scaling [9-12] — with the help of e.g. 
the NOON states [13]. 

Still, in all realistic experimental setups, deco- 
herence typically makes the relevant quantum fea- 
tures such as squeezing or entanglement die out very 
quickly. Recently, it has been rigorously shown for 




Figure 1. Quantum metrology with matrix product 
states. N parallel quantum channels act on the input 
state I*) inscribing parameter value 4> as well as causing 
local decoherence on each of the probes. Measurement 
IIj, on the output state allows to make an estimate 
(p{x) based on the measurement result x. The input state 
is a matrix product state obtained by the action of maps 
Aa^ on pairs of virtual D dimensional systems each pre- 
pared in a maximally entangled state \ipd)- 

optical interferometry with loss [14, 15], as well as 
for more general decoherence models [16, 17], that 
if decoherence acts independently on each of the 
probes one can get at best c/\/N asymptotic scal- 
ing of error — precision that is better than classical 
one only by a constant factor c which depends on 
the type of decoherence and its strength. 

This implies that in the limit of large number of 
probes, almost optimal performance can be achieved 
by dividing the probes into independent groups 
where only the probes from a given group are en- 
tangled among each other. Clearly, the size of the 
group that is needed to approach the fundamental 
c/y/N bound up to a given accuracy will depend on 
the strength of decoherence. Nevertheless, irrespect- 
ively of how small the decoherence strength is, for 
N large enough the; size of the group will saturate at 
some point and therefore asymptotically the optimal 
state may be regarded as only locally correlated. 

A natural class of states efficiently representing 
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locally correlated states are the Matrix Product 
States (MPS) [18-22], which have proved to be 
highly successful in simulating low-energy states of 
complex spin systems. Until now no attempt has 
been made, however, to employ MPS for quantum 
metrology purposes. Establishing this connection is 
the essence of the present paper. 

Basic quantum metrology scheme is depicted in 
Fig. 1. A'' probe input state l^*) travels through 
A'' parallel noisy channels which action is para- 
meterized by an unknown value (j). A measure- 
ment Tix is performed on the output density matrix 

= K®^ {\^){'^\) yielding a result x with probabil- 
ity p{x\4>) = Tr(/50na;). The estimation procedure is 
completed by specifying an estimator function 4'{x). 
Eventually wc are left with the estimated value of 
the parameter, </>, which in general will be different 
from (j). We denote the average uncertainty of es- 
timation by A(j!) = ((^ — 4>Y). where the average 
is performed over different measurement results x. 
The main goal of theoretical quantum metrology is 
to find strategies that minimize A(^. For this purpose 
one has to find the optimal estimator, measurement 
and input state. This in general is a difficult task. 

To simplify the problem one may resort to the 
quantum Cramer-Rao inequality [23-26] 

that bounds the precision of any unbiased estimation 
strategy based on k independent repetitions of an ex- 
periment. F{f)^) is the Quantum Fisher Information 
(QFI) written in terms of the so called symmet- 
ric logarithmic derivative (SLD) defined implicitly 
as: 2^^ = L^p^ + P(f>L(t>- For pure states the for- 
mula for QFI simplifies to F{\^^)) = 4{{i'^\i'^) - 
where j^-^) = The bound is 

known to be saturable in the asymptotic limit of 
A; ^ 00 in the sense that there exist a measurement 
and an estimator that yields equality in (1). The 
main benefit of using QFI is that since it does not 
depend neither on the measurement nor on the es- 
timator, the only remaining optimization problem is 
the maximization of F{p^) over input states. 

Since the optimal states in the regime of large 
number of probes N (not k\) may be regarded as 
consisting of independent groups, the Cramer-Rao 
bound may be saturated even for A; = 1 provided N 
is large enough [25, 27]. This makes the QFI an even 
more appealing quantity than in the decoherence- 



free case were some controversies arise on the prac- 
tical use of the strategies based on the optimization 
of the QFI [28, 29]. 

Maximization of QFI over the most general input 
states for large N may still be challenging and even 
if successful might not provide an insight into the 
structure of the optimal states. This is the place 
were MPS come in useful. 

Here for concreteness we focus on one of the most 
thoroughly analyzed and relevant model in quantum 
metrology namely the lossy interferometer. In or- 
der to treat uniformly both atom and photon in- 
terferometry, we will not specify the nature of the 
physical systems but will rather refer to abstract 
two-level probes, with orthogonal sates |0), |1) (e.g. 
a photon traveling in the lower/upper arm of the 
interferometer). Consider N distinguishable probes 
(e.g. photons prepared in different time bins), = 

Slii2...iw=o"H»2...iivKi«2 • • -iiv)- The parameter to 
be estimated is the relative phase delay </> a probe 
experiences being in |1) vs. |0) state. The decoher- 
ence mechanism amounts to a loss of probes where 
each of the probes is lost independently of the others 
with probability 1 — rj. As such, this is an example 
of a general scheme depicted in Fig. 1. Since the 
distinguishability of probes offers no advantage for 
phase estimation [30] it is more convenient to work 
with indistinguishable probes and move to the oc- 
cupation number basis where the general N probe 
state reads |^) = X^^o^nl^'-^"'^)' \n,N — n) 
represents n and N — n probes in states |0) and |1) 
respectively. Formally speaking, switching from the 
distinguishable to the indistinguishable probes case 
is equivalent to restricting the full N probe space to 
its fully symmetric subspace. The output state p^ 
can now be written explicitly as: 

N N-lo 

pt-=ii Ef^wj*r)(*n (2) 

;o=o ii=o 

where 

N-h 

(3) 

and pigi^ is a normalization factor which can be inter- 
preted as a probability to lose Iq, li probes in states 
|0) and ]1) respectively. 
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As the output state is mixed, QFI is not expli- 
citly given in terms of the input state parameters as 
it requires calculation of the SLD. As was shown in 
[31], an excellent approximation to QFI is obtained 
in the form of a weighted sum of QFIs for the pure 
states entering the mixture (2): 

N N 

F{h)^T.llPioMm't))- (4) 

/o=0 ii=0 

The approximation amounts to a statement that 
knowing how many probes were lost, additional 
knowledge of how many of them were lost being 
in the state |0) and how many being in the state 
1 1) does not improve the estimation precision appre- 
ciably. This simplifies the calculations significantly 
since QFI for pure states has an explicit formula — 
in our case F{\^f^)) = 4((*^"'i In^j^'^o'i) _ 
|(\IiWi|^|xj;Wi^|2-) .^ji-jj being the the excitation 
number operator n|n, N — n) = n|n, N — n). 

Direct optimization of formula (4) over the input 
state parameters q:„ involves A'^ + 1 variables. This 
approach was taken in [30, 31]. Here we want to 
show that by restricting ourselves to a class of low 
rank MPS we can achieve almost identical precision 
with optimization over very few parameters. 

Translationaly invariant MPS of N qubits with 
periodic boundary conditions are defined as 

|*)mps = ^ V Tr(^i, ...Ai^)|ii...ijv), 

VN ■ ■ n 
Zi...7Ar=0 

(5) 

where two matrices Aj are square complex matrices 
of dimension D x D where D is called the bond di- 
mension and Af is the normalization factor, Af = 
El..i.=o ® . . . ® A,)]. In op- 

erational terms, a MPS is generated by assuming 
that each qubit is substituted by a pair of D dimen- 
sional virtual systems. Systems corresponding to 
neighbouring particles are prepared in maximally en- 
tangled states \(Pd) = SjLi lijj) (Fig- 1) and maps 
Ai = X^^^^^i Ai^ct,fi\i){ct, 1^1 are applied to each pair 
of the virtual systems yielding eventually the MPS 
[21]. For the problem considered, the advantage of 
using MPS stems from the fact that for the bond 
dimensions D needed to approximate a given state 
using a MPS will typically be much smaller than A''. 
For the purposes of lossy phase estimation we can 
restrict ourselves to fully symmetric states for which 
a sufficient condition is that the matrices Aq, Ai be 
diagonal. Even though for a given D a larger class 



of symmetric states could be generated with more 
general matrices, we have checked that this does not 
provide any computational benefit in our case. As a 
results the MPS are parameterized with 2D complex 
numbers instead of A' -|- 1 needed to specify a gen- 
eral symmetric state of A^ qubits. In the standard 
occupation number basis the MPS has the form: 

|*)mps = ^f2 ^/(f) Tr(A^Af-")|n, A^ - n). 

(6) 

Whenever possible one should compute the quant- 
ities of interest directly on the Ai matrices. For- 
tunately, thanks to the simple form of Eq. (4) it is 
indeed possible and there is no need to go back to 
the standard state description as would be the case 
with the exact formula (1). 

Fig. 2a illustrates the precision obtained using 
MPS for the case of relatively small losses rj = 0.9. 
As one can see, the AIPS approximation is excel- 
lent. In particular, the upper-right inset shows that 
already D = 5 is sufficient to obtain less than 1% dis- 
crepancy for A^ < 500. We have confirmed this ob- 
servation for different rj and observed that for higher 
losses (lower 77) lower D are required to obtain a 
given level of approximation for a particular A^ — an 
effect that should be much more spectacular for lar- 
ger A^ reflecting the fact that stronger decoherence 
diminishes the role of quantum correlations. 

Moreover we have observed that optimal matrices 
Ai have the same diagonal values which 
are ordered complementarily — the highest in Aq is 
paired with lowest one in Ai etc. The higher is A' 
the closer the diagonal values approach each other 
as can be seen from the lower- left inset on Fig. (2). 
This confirms the intuition that with increasing N, 
the optimal states are becoming less distinct from 
the product state — all diagonal values of Ai equal. 

The peciiliarity of phase estimation is that in the 
decoherence-free case optimal QFI is achieved for the 
NOON state which, even though has non-local correl- 
ations, is an example of an MPS with D = 2. This 
makes the MPS capable of approximating the op- 
timal states very well even for low loss and small A" 
[N < 1/(1 — 77)] ^'^n ability that in general will not 
hold for other estimation problems. 

Until now we have focused our analysis only on 
the QFI and avoided discussion of any particular 
measurement scheme. Taking less fundamental but 
more operational approach, one may consider a con- 
crete measurement scheme were a particular observ- 
able O is being measured. Simple error-propagation 
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Figure 2. (color online) Log-Log plots of the phase estimation precision with losses {rj = 0.9) as a function of number 
of probes A'^ optimized over input states in (a) the QFI approach were Atp = with F given by Eq. (4) and 
(b) the Ramsey interferometry with A<j> given by Eq. (7). The color curves correspond to matrix product states of 
different bond dimensions D: D = 1 (green, no correlations: A(j} = l/^rjN), D — 2 (purple), D — 3 (red), D — 4 
(orange), D = 5 (yellow). The blue curves correspond to the optimal states obtained with brute force optimization, 
with extrapolation for higher A'^ (dashed) obtained using scaling formulas of [14] (a) and optimization over the spin- 
squeezed states [33] (b). The black line represents the asymptotic bound A(f> — ^ [14-17]. The upper-right 

insets depict the bond dimension D needed to reach the optimal precision with at most 1% discrepancy: r) — 0.9 
(solid), rj — 0.3 (dashed). Lower- left insets present normalized absolute values of diagonal matrix Ai elements for 
bond dimension D = 4 and rj — 0.9. 



formula for (f> yields A(/) = AO/\'^^^\. In the Ram- 
sey spectroscopy setup [32], or equivalently in the 
Mach-Zehnder interferometer with photon number 
difference measurement, one effectively measures a 
component of the total angular momentum operator 
J of N spins 1/2 — if a qubit |0), |1) is treated as a 
spin 1/2 particle [11]. If the phase dependent rota- 
tion U = e^'^'^' is being sensed by the measurement of 
the Jx observable, the explicit formula for estimation 
uncertainty at the optimal operation point <j) ~ cal- 
culated for from Eq. (2) reads 



A0 




(7) 



Search for the optimal state amounts to minimiz- 
ing the above quantity. Since it depends only on 
first and second moments of J it is simple to imple- 
ment numerically using MPS. Results are presented 
in Fig. 2b. 

It is clear that MPS are capable to capture the es- 
sential feature of the optimal states — the squeezing 
of the Jx — with relatively low bond dimensions D. 
Moreover, the upper-right inset indicates that the re- 
quired bond dimension D is reduced much more sig- 
nificantly with increasing decoherence strength than 



in the QFI approach. The lower-left inset confirms 
again that the structure of the optimal states gets 
closer to the product state structure with increasing 
N. 

In summary, we have shown that MPS are very 
well suited for achieving the optimal performance 
in realistic quantum mctrological setups and may 
reduce the numerical effort while searching for the 
optimal estimation strategies. Even though we 
have based our presentation on a single model of 
lossy phase estimation we predict this conclusions 
to be valid in all mctrological setups were deco- 
herence makes the asymptotic Heisenberg scaling 
unachievable — the intuitive argument being that no 
large scale strong correlations are needed to reach the 
optimal performance. An intriguing open question 
remains: is it possible, as it is in many-body physics 
problems, to obtain an exponential reduction in nu- 
merical complexity thanks to the use of MPS. This 
is not possible when the optimal states are known 
to be symmetric, as in the lossy phase estimation, 
since then the Hilbert space dimension grows only 
linearly with N. In problems, however, were distin- 
guishability of probes is essential as e.g. in reference 
frame estimation [34], MPS might demonstrate their 
full potential. 
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